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d ■ Abstract 

O^i We discuss the impact of a finite effective range on three-body systems interacting through a 



large two-body scattering length. By employing a perturbative analysis in an effective field theory 
well suited to this scale hierarchy we find that an additional three-body parameter is required for 
consistent renormalization once range corrections are considered. This allows us to extend previ- 
ously discussed universal relations between different observables in the recombination of cold atoms 



\^J I to account for the presence of a finite effective range. We show that such range corrections allow us 



to simultaneously describe the positive and negative scattering-length loss features observed in re- 
combination with "^Li atoms by the Bar-Ilan group. They do not, however, significantly reduce the 



"iJ I disagreement between the universal relations and the data of the Rice group on ^Li recombination 



at positive and negative scattering lengths. 
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I. INTRODUCTION 

In the 1970s Vitaly Efimov showed that the non-relativistic three-body system with two- 
body scattering length, a, much larger than the range of the underlying interaction displays 
universal properties that are independent of the details of the interaction [l|,l2|. In particular, 
in the unitary limit \a\ — )■ oo, there exists a tower of three-body states {trimers) with a 
geometric spectrum, 

4") = (e-2"/^«)"-"*n2«2/m, (1) 

where k* is the binding wavenumber of the trimer labeled by n*, m is the mass of the 
particles, and sq = 1.00624. Eq. ([1]) and other universal predictions have received significant 
attention recently in both atomic physics and nuclear physics. In the latter case, the three- 
nucleon system and certain halo nuclei display features reflecting their status as (effective) 
three-body systems which are near the unitary limit (for reviews see |3|, |j]). 

In atomic physics, the three-body recombination rate of ultracold atoms is an observable 
that displays the discrete-scale invariance of Efimov systems: it is proportional to a^ times 
a function that is log-periodic in a [3]. Several recent experiments have demonstrated the 
existence of such Efimov physics by measuring the recombination rate as a function of the 
scattering length. For example, Ref. p] measured the three-body recombination rate for 
^Li atoms. Gross et al. found a recombination minimum at a scattering length a = a*o ~ 
llGOa^, as well as a rate enhancement when a = a'^ ^ —264(11)0^. (Here and below qb 
denotes the Bohr radius.) 

The ratio a'^/a^o is very close to the "universal" result |3|, |6|: 

■(2n + l)7r 



— exp 



a*o " — > — 0.21a^,o. (2) 
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However, the interactions between ^Li atoms have a natural, van der Waals, length scale 
~ lOOa^, so significant corrections to this universal prediction for a'^ are expected. In par- 
ticular, the impact of a finite effective range on the scattering-length dependence of Efimov 
physics has been discussed in two recent papers [71, 181]. In contrast to those works, here we 
include the effect of a finite effective range perturbatively and analyze the consequences for 
universal predictions in the three-body sector. We do this by setting up a perturbation the- 
ory around the unitary limit, and thus organize the corrections that manifest the differences 



between different large-scattering-length systems. These are effects beyond universality, and 
their inclusion allows us to extend the reach of Efimov's ideas. The small parameters in this 
effective field theory (EFT) expansion are ^/\a\ and ik, where k denotes the momentum 
scale of the problem under consideration, and is ~ l/|ct| for the recombination processes 
that are our concern here. This "short-range EFT" applies to all non-relativist ic systems in 
which |a| is much larger than the range of the underlying interaction i. Part of the dynamics 
at scale d, enters observables via the two-body effective range, r^, and can be straightfor- 
wardly accounted for. But, other effects due to the finite range of the inter- atomic force 
get encoded in new three-body EFT parameters. These must be included in the calculation 
if it is to contain all effects up to a given order in the EFT. In this way we systematically 
approximate the dynamics of any finite-range interaction that gives a large scattering length 
a. The main result of our EFT analysis at next-to-leading order (NLO) is that, if (and only 
if) scattering-length-dependent observables are considered, one such additional three-body 
parameter must be included in the calculation in order to guarantee the accuracy of the 
EFT's predictions. To demonstrate the implications and limitations of this result, we use 
this approach to calculate the three-body recombination rate of ^Li in the hyperfine states 
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relevant to the experiments performed by Pollack et al. [9| and Gross et al. [5|. 

II. EFFECTIVE FIELD THEORY 

EFTs are a standard tool for calculating low-energy observables in systems with a sep- 
aration of scales. Here we will consider EFT for particles interacting solely via short-range 
interactions. Another successful example of an EFT is chiral perturbation theory (xPT) 



that identifies pions as the Goldstone bosons of low-energy QCD (see Ref. lOj for a recent 
review). At the heart of every EFT is a Lagrangian, which contains all possible operators 
allowed by the underlying symmetries. The short-range EFT applies to non-relativistic par- 
ticles with a large scattering length and contains only contact interactions. In this EFT we 
have, at leading order in an rs/\a\ (or, equivalently for our purposes, ^/|a|) expansion 

£ = V^t U -f!\^- ^{^HY - ^{^W + ■■■, (3) 

with ip our matter fields, and Co a two-body parameter fixed by the physical scale a. At 
leading order the two-body amplitude is obtained by summing all two-body diagrams that 



contain only the four-boson operator proportional to Cq. Requiring that the two-body 
amplitude has the physical scattering length leads to a renormalization condition for Cq: 

1 m mK , , 

if cutoff regularization with a cutoff A is employed. 

To go beyond leading order an ordering scheme derived from the requirement that each 
new order of the calculation gives contributions to observables that scale with powers of 

nn 

rs/|a| or krg is required |11| . |12| . This extends Eq. ([3]) to a predictive framework in which 
a well-defined class of terms has to be evaluated at every order in this small-parameter 
expansion. 



III. THE THREE-BODY SYSTEM 

In Refs. [l3|, llj] Bedaque et al. showed that this EFT is well-suited to describing the 
three-body problem near the unitary limit, and facilitates derivation of Efimov's results. In 
particular, they showed that a three-body parameter Dq must be present at LO in order to 
obtain renormalized quantities. 

This result was obtained by using the Lagrangian above to derive an integral equation for 
atom-dimer scattering. The S'-wave projected amplitude, ^o, for scattering from relative- 
momentum state k into relative-momentum state p, at energy E, is given by 



Aoip.k^E) 
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with e a positive infinitesimal. This equation (with Hq = 0) is known as the Skorniakov- 
Ter-Martirosian (STM) equation [l5'|. The additional term ~ Ho(A) can be thought of 
as a three-body force. It is proportional to Dq/Cq, and in general is a function of the 
integral-equation cutoff. The value of Hq is fixed by reproducing the value of one three- 
body observable, e.g. the binding energy of a particular three-body bound state at some 
fixed scattering length. This relates Hq to a physical scale in the three-body system. A 
common choice for this scale is k*, which can be thought of as the binding momentum of 
one three-body bound state in the limit \a\ — ?■ oo. Once k* — or equivalently Hq — is fixed. 



LO predictions follow for all other three-body observables, including observables at other 
scattering lengths which satisfy \a\ ^ L 

Yoi p\a\,p/k ^ 1, the asymptotic form of the amplitude Ao{p, k) is ~ p^^^o-i^ Corrections 
are suppressed by powers of l/(|a|p) or £p |l6|: 



Ao{p,k) =N'{k/k)\ -sin 
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Here M{k/h) denotes a normalization factor and C_i is a complex number that can be 



obtained from the Mel 
equation in Eq. ( 



in transform of the asymptotic form of the kernel of the integral 
161]. The need for a three-body contact interaction at leading order 
implies that — at fixed a — all observables in the three-body sector can be described by one- 
parameter correlations. For example, at fixed scattering length, particle-dimer scattering is 
correlated with the three-body binding energy. Mapping out such a correlation corresponds 
to varying K^,, or Hq^ keeping Cq fixed. One can, however, also vary the two-body contact 
interaction and thus change a, keeping the three-body counterterm fixed. This is what is 
expected to happen in atomic systems close to a Feshbach resonance. There, a small variation 
in the magnetic field changes the scattering length dramatically, but is not expected to affect 
the short-distance physics that determines the value of Hq. 

Since Hq (and hence k*) does not change as a is varied we have straightforward relations 
between the scattering lengths at which certain features in atomic recombination are ob- 
served. For example, the scattering length a* at which the trimer binding energy crosses 
the atom-dimer threshold is related to the binding momentum of the corresponding trimer 
at unitarity by a^, = 0.071/k* [3|]. 

Such relations can be listed for all relevant quantities and then lead to results such as 
Eq. ([2]). These equations will, however, be modified in presence of a finite effective range. 
In the next section we consider these effects. 



IV. BEYOND LEADING ORDER 



To include higher-order corrections we follow and extend^ Hammer & Mehen's analy- 
sis [ITI] and include the operator associated with the two-body effective range, r^, perturba- 



^ In Ref. 171 one term stemming from two-body scattering was neglected in the analysis. It will be taken 
into account here. 



tively in the calculation of observables. This treatment is valid provided that \rs\ ^ i '^ a. 
To calculate effects ~ r^ we consider the piece of the two-particle propagator that contains 
one insertion of effects beyond LO: 

^i(|q|; go) = - . . (7) 

^ — l/a + y —mqQ + q^/4 — le 

Diagram (a) of Fig. [1] then corresponds to a single insertion of the propagator (jTj) in the 
three-body system. It was shown in Ref. [l7| that the corresponding integral displays a 
divergence that has to be canceled with an energy- independent three-body force ifi(A). 
Diagrams (b)-(e) in Fig. 1 have to be evaluated for this purpose. 
When this is done we obtain: 

A,{k,k) = ':iAo{k,k)+ [ dqq'-^V,iq;E)Alik,q) 
aHi{A) 



7r3A2 



4vr2 , f^. 2 Ao{k,q) 

' ' dqq 



(8) 



^0 -l/a+ ^3q^/A-mE-ie 

where the first term comes from the NLO renormalization of the dimer field. Inserting 
Eq. (ini) into the integral in the first line of this expression makes it clear that the term ~ Hi 
on the second line has to absorb two independent divergence structures, one of which is 
proportional to l/a. Therefore we must write 

i/i(A) = i7io(A) + l^niA) ■ (9) 

Considering only fixed scattering length does not introduce an additional three-body param- 
eter, since the total Hi{A) is fixed by the renormalization condition that the NLO amplitude 
has zero effect on the observable used for renormalization at leading order. However, if we 
desire predictions at NLO for observables as a function of a we need to know the relative 
sizes of if 11 and Hiq. 

The existence of such an a-dependent three-body force can be understood by considering 
the Lagrangian shown in Eq. ([3]). This Lagrangian is constructed by writing down all terms 
allowed by the symmetries of the system. This includes terms proportional to powers of the 
small scales inherent to the problem: in our case, l/a. Beyond LO the three-body term in 
Eq. ([3]) must be augmented to: 

Do(^V)' + — (^V)' + ^{i^W + ■■■■ (10) 

a a^ 
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FIG. 1: Diagrams for the NLO amplitude of atom-dimer scattering. The crossed double line 
denotes an insertion of the range propagator, ^i(|q|;(?o); the shaded blob denotes an insertion of 
the LO three-body amplitude, ^O; and the square denotes an insertion of the NLO three-body 
force, Hi. 

This exemplifies the benefit of any EFT: all operator structures allowed by the under- 
lying symmetries will contribute at an order in the small-parameter expansion that can be 
determined by careful analysis. It furthermore relates experiments with ultracold atoms to 
a problem encountered in the application of xPT to nuclear physics. In this description of 
low-energy nuclear processes using only baryonic and pionic degrees of freedom, the small 
scale analogous to 1/a is the pion mass m.^, and operators proportional to powers of the pion 
mass are a mandatory part of expressions for low-energy nuclear observables. For example, 
at 0{m1), the nucleon-nucleon (NN) potential calculated in xPT can be written as 

V = [Cs + Dsml] + [Ct + DTml]ai ■ a^ + Vope + V^% , (11) 

(2) 

where Vqpe is one-pion exchange and V^p^; is the "leading" two-pion exchange contribution 



to the A^A^ interaction. Both Cs/t and Ds/t are required for consistent renormalization |18 |. 
The relative values of these two are not required for the calculation of observables measured 
in the laboratory, i.e. at fixed m^. However, as lattice QCD attempts to predict and explain 
few-hadron properties, the determination of pion-mass dependent coefficients becomes a 
crucial feature of future progress [l9|. 



V. THREE-BODY RECOMBINATION 



A. The Lithium-7 System 



The new three-body counterterm Hn (= Dn) is only relevant if systems with a variable 
scattering length are considered. Experiments with ultracold atoms are therefore the ideal 
place to explore its impact. Here we focus on ^Li, where we have information on how both 
the scattering length and the effective range vary with magnetic field, B. The dependence 
of a and r^ on the magnetic field as presented in Ref . J5| is shown in Fig. [2] for two different 
hyperfine states [20j. 

We will use our EFT, with this two-body input, to discuss two recent experiments which 
measured three-body recombination rates in systems of ^Li atoms. First, Gross et al. have 
measured the three-body recombination of atoms in the \F = 1 mp = 0) hyperfine state |5|. 
They found a recombination minimum at a > and a recombination maximum at a < 
whose relative positions are quite well described by universal (r^ = 0) predictions. Second, 
Pollack et al. have measured the recombination rate of atoms in the \F = 1 mp = 1) state 
and found several recombination features associated with few-body universality |9(] . In these 
data there was a systematic deviation by a factor of two from the universal prediction in 
the ratios of features on the positive and negative scattering-length sides of the Feshbach 
resonance. 

Using the diagrammatics discussed above we perform an EFT calculation of the three- 
body recombination rate into dimers with binding energy ~ l/ma^, for positive scattering 
length. Since our calculation does not include dimers bound by ~ l/mi"^, we cannot compute 
the recombination rate for negative scattering length. We can, however, determine the 
position of recombination maxima at a < 0, by calculating the scattering lengths for which 
the binding energy of a trimer becomes zero. 

The recombination length, ps, is obtained from the LO and NLO elastic atom-dimer 
scattering amplitudes, ^o and Ai, via (in units where h = 1): 



P3 




{Ao + Ai)(o,^;o\ 



(12) 



with 7 = y/mB2, and B2 the binding energy of the atom-atom dimer (7 = 1/a at LO). Figure 
[3] shows our results for p^ as a function of the scattering length a for the \F = 1 mp = 0) 
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FIG. 2: (Color online) Magnetic-field dependence of the scattering length and effective range in the 
hyperfine states relevant to the experiments by Gross et al. (left panel) and Pollack et al. (right 
panel). The solid lines denote the scattering length and the dashed lines denote the effective range. 



hyperfine state. The dashed line denotes the LO result renormalized to the recombination 
maximum a'^ on the a < side. This LO calculation predicts a*o > 12000^, c.f. the 
measured a*o ~ 11600^. Because, as discussed above, consistent renormalization at NLO 
requires us to choose two different three-body observables to fix the counterterms Hqq + Hiq 
and Hii, we can describe both the position of this recombination minimum and that of the 
observed maximum at a < 0. The solid line displays that NLO result, renormalized to a'^ and 
the recombination minimum a*o determined by Gross et al. [5]. The shaded area denotes 
the region where |rs|/a > 0.5 and convergence of the EFT expansion is expected to be 
slow. The squares give the experimental data with the corresponding errors. We emphasize 
that we did not strive for a detailed reproduction of the experimental data here, since our 
predictions account for neither the effects of deep dimers nor those of finite temperature. The 
inclusion of such physics improves the overall agreement between experiment and theory. 
To demonstrate this we also present the LO result with deep-dimer effects included: it is 
the dot-dashed line in Fig. [3l This shows that the positions of the loss features which we 
are focusing on in this letter are not affected by recombination into deep dimers. They are 
also not affected by finite temperature |2l| . 

It is gratifying that we can explain Ref. |5|'s measurements of both a' and a^,o and still 
obtain results consistent with the EFT expansion, but neither of these is an NLO prediction 
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FIG. 3: (Color online) The three-body recombination length of ^Li in the \F = 1 nip = 0) 
hyperfine state as a function of the scattering length. The dashed line denotes the LO result where 
the three-body parameter is determined by the recombination maximum that has been measured 
experimentally in Ref. [5[. The dot-dashed curve is the LO result including deep-dimer effects with 
fitting parameters as in Ref. [5|]- The solid line denotes our NLO result (which does not include 
deep dimers) 



of the EFT. We now use our NLO calculation to predict the scattering length at which the 
atom-dimer resonance occurs, a*, i.e. the atom-atom scattering length at which the atom- 
dimer scattering length diverges. The LO result for a^ is a* = — 1.03a'^, which puts a^: in the 
region where It's! /a ~ 0.3. Therefore effective-range corrections to this observable can be 
large. The result for a^, thus depends strongly on the LO amplitude that is employed in this 
calculation and so on the observable that was used at LO as three-body input. We obtain 

a* = (271 - 105 + . . .)aB (13) 

when a'^ is used at LO and a*o is used in addition at NLO. Alternatively, we find 

a, = (257 - 4 + . . .)aB (14) 

when a^o is used at LO and a'^ is used in addition at NLO. In either case, our results suggest 
that a* is shifted to smaller values than the LO prediction. The significant difference between 
the NLO prediction in the two different renormalization schemes provides an estimate of 
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FIG. 4: (Color online) The three-body recombmation length of '^Li in the |-F = 1 mp = 1) hyperfine 
state as a function of the scattering length. The dashed line is the LO result, with the three-body 
parameter determined by the recombination maximum on the a < side at a} = —63010^ . 
The dot-dashed curve is the LO result including deep-dimer effects with the deep-dimer fitting 
parameter of Ref. [9]. The solid line is our NLO result with the second three-body parameter 



adjusted to the second recombination maximum at a 



/(2) 



-298aB. 



effects at NNLO and beyond. In fact, we can expect that NNLO corrections are larger 
for the result (TT3|) . since the LO renormalization point is further from the unitary limit 
there. But, to encompass both ( 1131) and (JH]), as well as provide an uncertainty based on 
conservative estimates of NNLO effects, we quote: 



a* = (210±44)aB. 



(15) 



We have performed a similar analysis for the results obtained by the Rice group in 
Ref. |9| . In Fig. H] we show our results for the recombination length of ^Li atoms in the 
\F = 1 nip = 1) state as a function of the scattering length. The LO result (dashed line) is 
renormalized to the first recombination maximum a} = — GSOla^ on the a < side. The 
solid line shows the NLO result that is additionally renormalized to the second measured 

/(2) 

recombination maximum a* = —2980^. Having adjusted both three-body parameters, we 
can use them to calculate other observables such as 



K* 



(0.242 ± 0.030) • lO^^a 



B ' 



(16) 
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and the position of the atom-dimer resonance 

a, = (355.8±55.5)aB . (17) 

In the latter case a shift upward from the LO prediction ai = 295aB is seen. We can also 
predict the position of the second recombination minimum in Fig. H] to be: 

a,o = (1348 ± l51)aB • (18) 

The errors in Eqs. flTBl) - flTS]) were obtained by propagating the systematic and statistical 
errors quoted in Ref. [9|. The shift of the central value in fll8p from the LO result is only 26as, 
and so range effects ~ r^ and higher are certainly much smaller than the uncertainty due 
to the experimental input. The EFT result flTSj) therefore disagrees with the experimental 
result obtained in Ref. |9| 

a,o = (2676±67±128)aB , (19) 

by a factor of two, even after range corrections are included. (In Eq. (TT9!) the first error 
is statistical and the second due to a systematic uncertainty in the determination of the 
atom- atom scattering length.) Effects due to the finite effective range are therefore not 
responsible for the disagreement between the data of Ref. [9| pertaining to different sides of 
the Feshbach resonance and the predictions of universality. 

VI. CONCLUSION 

We have shown that at next-to-leading order in the EFT for systems with large scat- 
tering length an additional three-body input is required to describe physical observables 
as a function of the two-body scattering length. This parameter is needed for a complete 
NLO description of such observables if \rs\ ~ i. It encodes the effect of dynamics at dis- 
tances comparable to the range of the underlying two-body force on the scattering-length 
dependence of quantities in the three-body system. The necessity of this operator for renor- 
malization exemplifies the general tenet of all EFTs, that all operators which are allowed 
by the symmetries of the system will appear in the EFT Lagrangian. An expansion of the 
short- distance physics around the unitary limit (i.e. in powers of 1/a) makes clear that such 
an operator should be present at NLO in the short-range EFT. In this respect 1/a in this 
EFT plays a similar role to that of the pion mass in xPT. 
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Here we have analyzed the typical situation where the interaction range and the two- 
body effective range are of the same size: £ ~ \rs\. We expect the impact of our additional 
three-body counterterm to be smaller if £ ^ \rs\, i.e. the effective range is unnaturally large. 
This occurs, for example, close to a narrow Feshbach resonance. In that case the results of 
Ref. [Tj should hold. It is possible that the iteration of r^ to all orders can be systematically 
justified in such systems. In that case a three-body counterterm would not be required at 
LO of the EFT calculation [22;] and corrections to that "resummed LO" would scale as Vs/a 
and krg. 

We have shown that our approach can be applied to three-body recombination provided 
that the effective range is known. In the case of ^Li we calculated the shift in recombination 
features due to an effective range that is a function of the magnetic field. We found that 
deviations from universal predictions in the Bar-Ilan experiment can be explained by these 
effects but that they are not sufficient to explain the inconsistencies across the resonance in 
the experiment carried out by the Rice group. 



t 



The data of Ref. |9|] is not consistent with ~ r^ corrections to universality in spite of the 
presence of the extra three-body parameter at that order in our calculation. It has been 
suggested that this apparent violation of universality is a result of the conditions in the Rice 
experiment. For positive a the Rice group used a BEC, while for negative a they used a 
thermal gas. In contrast, the Bar-Ilan group recently repeated the Rice ^Li experiment in 
the \F = 1, rriF = 1) state using a thermal gas on both sides of the Feshbach resonance. 
They found that features across the resonance were related by universality [23|. 

Our results show that effects proportional to r^ correct the universal relations displayed, 
e.g. in Ref. [SJ] in a way that improves the agreement with data in the case of ^Li atoms |5|. 
It would be very interesting to apply this framework to data on three-body recombination of 
^^^Cs atoms. The existing data there, though, is rather different, since the features observed 
in experiments with ^^'^Cs at a > and a < are connected through the region where 



a = 24] • I^ the vicinity of this scattering-length zero Vg/a diverges, and higher-order 



corrections to the short-range EFT cannot be reliably calculated. The connection between 



a* 25|, a'^, and recombination minima at a > lOOOa^ can, however, be addressed within 
EFT. Given information on how the ^^^Cs two-body scattering length and effective range 
vary with magnetic field, we can make quantitative predictions for the impact of the range 
on three-body recombination features seen in experiments with cold gases of Cesium atoms. 
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